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SUMMARY 


resistivity  maps  from  which  sites  may  be  selected  lor  m J * C°nt°Urcd  int0 

(e.g.,  drilling).  morc  lnlenslve  ground  explorations 

intJJ^r^a^a  y^  °'  V which  ma*  aid  in  the  data 

variables  considered  were-  1)  contrast  in  resi  t surface  of  lhe  an°maly.  The  independent 
rounding  earth  2)  anomaly  szeni  V bCtWCC"  'hc  aild  lhe  «"■ 

maly,  and  (2)  a more  resistive  anomaly  Two  examnl«  fr„m 
mansions  ot  a ground  surface  anomaly. 
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NUMERICAL  STUDIES  TO  AID  INTERPRETATION 
OF  AN  AIRBORNE  VLF  RESISTIVITY  SURVEY 


Steven  A.  Arcone 


INTRODUCTION 

Background  and  objectives 

A knowledge  of  ground  resistivity  is  important  to 
many  engineering  projects,  since  it  can  provide  infor- 
mation for  determining  the  locations  of:  1 ) possible 
sites  for  providing  electrical  grounding  capability  (low- 
resistivity  sites),  and  2)  possible  sites  for  providing 
near -surface  deposits  of  sands,  gravels  and  bedrock 
suitable  for  construction  purposes  (usually  high- 
resistivity  sites). 

To  find  these  locations  rapidly  and  to  aid  in  charac- 
terizing the  ground  conditions  of  a region  rapidly,  an 
airborne  geophysical  system*  was  developed  (Barringer 
1972)  and  evaluated.  This  system  computes  an  appar- 
ent earth  resistivity  from  the  tilt  (from  vertical)  of  the 
electric  field  vector  of  low-frequency  (200-400  kHz) 
and  very -low-frequency  (10-30  kHz)  radiowaves.  Sur- 
veys are  usually  made  along  equispaced  flight  lines  and 
the  results  are  then  extrapolated  into  the  form  of  re- 
sistivity contour  maps.  These  maps  are  commonly 
superimposed  upon  a topographic  or  aerial  photomo- 
saic base  from  which  areas  of  interest  can  be  selected 
for  further  ground  tests. 

The  criteria  for  choosing  an  area  are  based  on: 

1 ) the  proximity  of  the  area  to  a relevant,  planned  con- 
struction site,  2)  the  contrast  between  the  resistivities 
of  an  anomalous  area  and  a more  common  area,  and 
3)  the  extent  of  a specified  resistivity  contour  desig- 
nated as  characteristic  of  the  material  or  conditions 
sought. 

The  importance  of  item  1 is  self-evident.  Items  2 
and  3 are  obviously  affected  by  the  resistivities  of 
individual  locations  and  the  regional  extent  of  the 
materials  present.  Item  3 is  specifically  a method  for 

* The  principles,  techniques  and  practice  of  this  system 
are  described  in  Hoekstra  et  at.  (1974)  and  will  be 
briefly  reviewed  in  the  next  section. 


estimating  the  actual  size  of  the  material  present.  How- 
ever, previous  field  data  (Hoekstra  et  al.  1975)  indi- 
cated that  a great  deal  of  smoothing  of  ground  resis- 
tivity information  occurs  even  at  minimal  flight 
altitudes  of  75  meters.  Therefore,  an  objective  was 
established  to  develop  procedures  for  better  esti- 
mating the  individual  effects  of  flight  altitude,  and  the 
areal  extent  and  relative  resistivity  of  an  anomalous 


scope 

A VLF  plane  wave  was  assumed  incident  upon  a flat 
earth.  A resistivity  contrast  was  established  over  a cir- 
cular portion  of  the  surface.  The  secondary  electric 
fields  that  resulted  from  this  change  were  then  calcu- 
lated at  a specified  survey  altitude  and  superimposed 
upon  the  general  plane  wave  field.  The  resulting  field 
ratios  above  the  anomaly  were  converted  to  apparent 
resistivity  values.  These  values  were  used  to  demon- 
strate the  result  of  the  ground-to-air  transformation 
and  to  determine  the  approximate  size  and  intensity  of 
a ground  anomaly  from  a knowledge  of  the  ground 
resistivity  values  obtained  from  airborne  surveys. 


PRINCIPLES  OF  VLF  RESISTIVITY 
SURVEYING 

This  section  briefly  reviews  the  principles  and  tech- 
niques involved  in  ground  resistivity  surveying  at  very 
low  frequency  (VLF).  A more  complete  discussion 
of  the  electromagnetic  theory  involved  may  be  found 
in  books  by  Wait  (1962)  and  Watt  (1967).  Surveying 
technique  and  results  are  covered  in  articles  by  Hoekstra 
et  al.  (1974, 1975),  Palacky  and  Jagodits  (1975), 
McNeill  et  al.  (1973),  Frischknecht  (1971)  and  Keller 
et  al.  (1970).  Discussions  of  the  resistivity  of  earth 
materials  are  presented  in  Ward  (1967),  Keller  and 
Frischknecht  (1966)  and  Parkhomenko  (1967). 





f igure  I.  I wo  radiowave  propagation 
modes  in  the  earth-ionosphere  wave- 
guide (see  text  for  explanation  ot  sym 
hols).  The  electromagnetic  field  c ompo- 
nents  are  characteristic  ot  radiation  from 
a vertical  electric  dipole  or  monopole. 


Electromagnetic  theory 

Ideally  suited  for  probing  the  resistivity  structure 
of  the  first  few  hundred  feet  of  the  earth’s  crust  are 
VLf  (10-30  kHz)  radiowaves.  Several  VLF  trans- 
mitters are  located  throughout  the  world  operating  at 
radiated  power  levels  between  0.1  and  I MW.  At  any 
point  on  the  globe,  at  least  one  ol  these  transmitters 
can  be  monitored  because  of  this  high  power  and  the 
confinement  of  the  radiation  to  the  earth-ionosphere 
waveguide. 

Within  this  waveguide,  radiation  propagates  in  two 
important  modes,  the  surface  wave  mode  and  the  skv 
wave  mode,  as  illustrated  in  I igure  1 . The  surface 
wave  mode  is  so  named  because  the  radio  wave  field 
strengths  attenuate  exponentially  with  height  above 
the  earth's  surface.  The  sky  wave  mode  is  the  result 
ol  many  "hops”  or  skips  oil  the  ionosphere  and  earth. 
At  the  earth's  surface  the  electromagnetic  vectors  ol 
both  these  waves  are  similarly  oriented.  Referring  to 
the  figure,  f represents  the  electric  field  components 
and  H the  magnetic  field  component  in  local  \,  v,/ 
right-hand  Cartesian  coordinates,  a and  ^ refer  to  the 
incident  and  grazing  angles  upon  the  ionosphere  and 
earth  respectively.  I or  sky  waves  this  field  configura- 
tion results  from  the  vectorial  addition  ot  the  incident 
and  reflected  waves;  it  is  referred  to  as  transverse 
magnetic  and  is  associated  with  the  vertical  electric 
dipole  and  monopole  types  ot  antennas  which  are 
commonly  used. 

Information  about  the  resistivity  ol  the  earth  may 
be  obtained  bv  measuring  the  complex  ratio  ol  / JL 
(termed  the  wavetilt  W)  or  of  L x///y  (termed  the  sur- 
face impedance  /,).  At  distances  greater  than  about 
8(X)  km,  i v -»  O''  for  VLI  waves  and  the  radiation  is 
essentially  planar.  The  wavetilt  at  the  surface  of  a 


homogeneous  earth  then  can  be  well  approximated 
as 

W = s/Jcac0p  ( | ) 

which  may  also  be  expressed  as 
W r sJu>eQp  exp  (/45°) 

where  2rr / where  t is  the  operating  frequency, 

in  Hertz 

c0  8.85x10"' 2 farads/meter 
p resistivity,  in  ohm-meters 

i V-T 

I or  all  simply  layered  cases,  /s  is  simply  related 
to  If  by  the  equation 


7s  v'M0/«o  W (2) 

where /Uq  4rrxl0~  henrys/metcr. 

Numerically  /s=  377  IV.  I quations  I and  2 assume 
that  displacement  currents  are  of  little  consequence 
at  VLF , which  is  reasonable  lor  most  earth  materials 
(lor  which  resistivities  are  commonly  below  10,000 
ohm-m). 

Resistivity  p is  then  obtained  by  simply  inverting 
eq  I or  2 and  is  calculated  using  the  formula 

P - |M’|2/co(0  (3) 

and/or 

p -|/s|2/ujp0 


2 


(4) 


which  gives  the  correct  value  of  ground  resistivity  for 
a homogeneous  earth. 

The  phase  angle  of  45°  given  in  eq  I is  diagnostic 
ol  homogeneity  to  at  least  the  depth  6 given  by  the 
formula 


6 is  commonly  referred  to  in  the  literature  as  the  shin 
depth,  the  depth  at  which  the  refracted  wave,  shown 
in  Figure  1 , attenuates  to  e-1  of  its  original  value  of 
the  surface.  Its  derivation  may  be  found  in  most 
textbooks  |e.g.,  Stratton  (1941)  and  jackson  (1962)] 
on  electromagnetic  theory.  At  VLF,  6 ranges  from 
about  36  m in  a 100-ohm-m  material  to  about  360  m 
in  a 10,000-ohm-m  material. 

When  resistivity  layering  begins  above  a depth 
equal  to  the  skin  depth  of  the  upper  layer,  the  magni- 
tude and  phase  of  IF  (or  /$)  changes  according  to  the 
resistivities  and  any  layer  dimensions  involved.  A 
measurement  of  the  ground  resistivity  is  still  converted 
from  either  It'  or  /s  as  in  eq  3 and  4,  but  only  gives 
some  representative  value  lot  the  composite  geoelectric 
structure.  Therefore,  p must  be  termed  an  "apparent” 
resistivity  and  is  customarily  symbolized  as  pa.  Phase, 
which  must  always  remain  between  0°  and  90°  for 
vertically  stratified  homogeneous  layers,  then  becomes 
diagnostic  of  the  layering  present.  In  general,  the 
following  rules  apply: 

1)  Phase  angles  45°  indicate  that  resistivity  is  de- 
creasing with  depth. 

2)  Phase  angles  < 45°  indicate  that  resistivity  is  in- 
creasing with  depth. 

formulas  are  available  to  account  for  these  effects 
(Wait  1962),  but  unless  more  than  one  frequency 
(transmitter)  is  available  for  measurement,  these 
formulas  arc  ol  little  help  in  resolving  subsurface  layer- 


Measurement  techniques 

In  the  airborne  system  a pair  ol  horizontal  and 
vertical  electric  dipoles  are  used  to  obtain  a measure 
ol  the  wavetilt  magnitude.  / , is  selected  as  the  phase 
reference  against  which  both  in-phase  and  quadrature 
phase  components  of  L may  be  compared  to  obtain 
the  total,  complex  ff.  However,  Might  instabilities  of 
the  aircraft  usually  cause  the  vertical  field  to  couple 
into  the  horizontal  antenna,  thereby  producing  dis- 
tortions in  the  in-phase  component  of  IF.  Consequent- 
ly, only  the  quadrature  component  is  available  for 
computing  ground  resistivity.  This  computation  then 
introduces  a discrepancy,  since  it  must  assume  an 


arbitrary  phase  angle  ol  45°,  appropriate  tor  only  the 
homogeneous  model. 

The  formula  for  conversion  from  wavetilt  to  ap- 
parent resistivity  is  then 

2 Quad2  (Ex/Ez) 


The  relationship  between  the  actual  magnitude  ol  the 
apparent  resistivity  and  that  computed  from  the 
quadrature  wavetilt  component  is 

Pa  . =2p  sin20  (7) 

aatrborne 

where  0 is  the  phase  angle.  Consequently,  there  is  a 

system  emphasis  on  resistivity  values  for  materials 

nearer  the  earth’s  surface.  This  is  because  p 
. i . . . . _o  , . “airborne 

increases  when  0 is  greater  than  45  and  decreases 

when  it  is  below  45°  in  accordance  with  phase  guide- 
lines (1 ) and  (2)  on  this  page. 

Ground-based  systems  are  able  to  measure  both 
phase  and  amplitude  of  Zs.  If  is  not  measured  on  the 
ground  because  vegetation  strongly  disturbs  Ez  at  the 
earth's  surface. 


THEORETICAL  PROCEDURES 
Assumptions 

The  basic  assumptions  employed  are:  1 ) the  sur- 
lacc  is  flat,  2)  the  VLF  waves  arc  incident  at  a finite 
but  small  grazing  angle  0,  3)  the  VLF  waves  are  trans- 
verse magnetic  and  plane,  and  4)  only  E x changes  at 
the  surface  of  the  anomaly. 

The  final  assumption  is  based  on  the  evidence  that 
where  topographic  relief  is  not  severe  and  no  highly 
conductive  zones  (ore  bodies)  arc  present,  Ez  and  //y 
remain  relatively  stable,  while  only  E x varies  suf- 
ficiently with  resistivity  changes  to  allow  a qualitative 
differentiation  of  material  type  based  on  the  computa- 
tion of  apparent  resistivity. 

The  particular  cases  considered  assumed  that  the 
homogeneous  earth  value  0 = 45°  does  not  change 
over  the  anomaly.  This  is  not  necessary  for  the  com- 
putational program  to  be  discussed  and  any  value  may 
be  considered  for  either  the  earth  or  the  anomaly. 

Tlie  anomaly  chosen  is  represented  at  the  earth's  sur- 
face as  a circular  area,  although  this  too  is  not  neces- 
sary. It  is  only  chosen  lor  its  simplicity.  A discon- 
tinuous change  in  the  strength  of  E x occurs  upon  en- 
countering the  anomaly.  Although  this  is  theoretic- 
ally impossible,  it  is  a reasonable  approximation, 


U 
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Figure  2.  Uniform  earth  of  resistivity  p,  containing  an  anomalous  area  of  resistivity  p2.  The 
electromagnetic  field  vectors  refer  to  the  incident  (at  an  angle  O-left)  and  reflected  (right)  waves. 
E'xa  is  the  perturbed  electric  field  on  the  anomaly.  Surface  coordinates  are  primed  and  spatial 
coordinates  are  unprimed.  Both  coordinate  systems  have  coincident  origins.  The  radius  of  the 
circular  area  is  rB. 


since  the  results  of  ground  surveys  discussed  in  Blom- 
quist  (1970)  and  Hoekstra  et  al.  (1975)  have  shown 
that  rapid  changes  in  apparent  resistivity  can  occur 
over  distances  of  only  tens  of  meters. 

Mathematical  formulations 

Far  from  the  transmitter,  the  electromagnetic  fields 
of  a VLF  sky  wave  can  be  represented  at  and  above 
the  earth’s  surface  as  the  sum  of  incident  and  reflected 
waves  as  follows: 

Hy  = H0  exp  (~ik0x  sin 0 ) [exp  (ik0z  cos 0) 


where  H0  = incident  magnetic  field  strength,  in  amp- 
turns/m 

0 = complement  of  ^ in  Figure  1 
R = reflection  coefficient  which  depends  upon 
0,  co,  and  the  electrical  properties  of  the 
earth 

?0  = VFo/eo 
Af0  = Wp0e0  . 

A time  dependency  of  exp  (/cot)  is  assumed.  Ex  andfz 
are  expressed  in  volts/meter  and  are  related  to  the 
magnetic  field  through  the  vector  relation 


+ R exp  (-iko?  cost?))  (8) 

Ez  - H0Z0  si nO  exp  (-/'Aq*  si nO)  [exp  (/'AgZ  cosfl) 
+ R exp  (-ikfjZ  cosfl)]  (9) 

and 


E~=t-^-  VxH 


where  V is  the  gradient  operator  and  is 


V=?±  + i±  + k*- 
dx  'by  dz 


(11) 


(12) 


Fx  = HqZq  cos 0 exp  (-ikQx  sin 0)  [exp  [ik^z  cosfl) 
-R  exp  (ikfyZ  cosfl)]  (10) 


* These  equations  satisfy  the  homogeneous  (source- 
free)  Helmholtz  wave  equations,  either  V^H  + 

= 0 orV^E  + k\E  = 0. 


A ^ A . 

where  /,  /,  and  k are  unit  vectors  in  the  x,  y,  z directions 
respectively  for  a right-hand  coordinate  system.  This 
system,  the  unit  vectors,  and  the  incident  and  reflected 
fields  are  illustrated  in  Figure  2.  The  vector  A’0  in  the 
figure  refers  to  the  wave  propagation  direction  and  the 
unit  vector  n is  the  surface  normal. 

The  effect  of  a vertical  stratification  of  any  number 
of  homogeneous  layers  of  different  thicknesses  and 
resistivities  is  expressed  through  the  reflection  coefficient 
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R.  At  given  values  of  cu  and  6,  R may  be  expressed 
using  the  surface  impedance  Zs  as  follows: 

* = (!-«/(  l+fl)  03) 


Referring  again  to  Figure  2,  the  secondary^clectric 
fields  induced  by  a specified  distribution  of  E ' and 
H'  over  a finite  surface  area  are  found  from  the  vector 
integral  equation  (Stratton  1941,  p.  476) 


where 


0 = ZJZ0  cosfl. 


(14) 


f = - jjj - f J [ioJH0(nxH')±  + (z?x£') 


This  allows  /?  to  be  computed  from  any  complex 
value  of  Z%  that  may  be  determined  experimentally. 

Equations  8,  9,  and  10  now  represent  an  average, 
background,  plane  wave  field  upon  which  perturba- 
tions resulting  from  changes  in  Ex  are  to  be  super- 
imposed. Technically,  a spatial  change  in  only  £x, 
without  accompanying  changes  in  the  other  Held  com- 
ponents, violates  the  divergence  equation  V-E  = 0 for 
a source  free  electromagnetic  field.  However,  this  is 
assumed  to  be  a reasonable  approximation  based  on 
the  experimental  evidence  referred  to  in  the  previous 
section.  Therefore,  these  local  changes  arc  represented 
as 


C = (ZJZsK  05) 

where  prime  = a quantity  evaluated  at  the  ground 
surface 
z = 0 

Z a = new  value  of  surface  impedance 
£'xa=  new  surface  value  of  Ex  on  the  anoma- 
lous area. 

The  net  change  from  the  plane  wave  value  in  electric 
field  intensity  E"xa  on  the  surface  of  the  anomaly  is 
then 


The  secondary  fields  induced  by  this  local  change 
in  £'x  arc  usually  referred  to  as  ‘'scattered"  fields.  The 
classical  approach  to  scattering  problems  involves  find- 
ing an  equivalent  set  of  surface  currents  and  surface 
charges  to  replace  the  scattering  obstacle  in  question 
(Harrington  1961).  The  italicized  terms  are  nomen- 
clatures for  the  vector  and  scalar  products  of  the  sur- 
face normal  (n  in  Fig.  2)  in  question  with  the  associated 
surface  electric  and  magnetic  fields.  Since  the  surface 
electromagnetic  representation  of  the  resistive  anomaly 
has  already  been  assumed,  cq  8,  9,  1 0,  and  1 6 now  de- 
fine these  surface  fields.  For  an  idealized  flat  earth,  the 
surface  normal  coincides  with  the  z axis. 


x V*  + (n.E')V$]do'  (17) 

where 

4^  = exp  (~ik0r)  /r  (18) 

r = \/(x  - x')2  + (y -y')2  + z2  (19) 

da'  = an  elemental  area  of  integration 
x'  andy'=  surface  coordinates  for  the  area  of  integra- 
tion in  the  plane  defined  by  z - z = 0. 

The  spatial  coordinates  of  the  point  at  which  the 
fields  are  to  be  evaluated  are  defined  as  x,  y,  and  z. 
Equation  1 7 is  often  referred  to  as  the  vector  equiva- 
lent of  Kirchoff’s  scalar  diffraction  integral. 

Since  the  only  perturbation  considered  is  £xa,  the 
total  airborne  electric  field  is  then  the  sum  of  the 
incident,  reflected  and  scattered  anomaly  fields. 
Mathematically,  this  total  field  is 

E = ?EX  +kEz-± - j f (nxt"J  x V±da'  (20) 

x'  y' 


where  Ex  and  Ez  are  specified  by  eq  9 and  10.  Carry- 
ing out  the  differentiation,  it  is  found  that 


(21) 

where  /0  is  the  specified  altitude  of  surveying.  The 
quadrature  value  of  the  wavetilt  is  then  computed  from 
cq  21  and  converted  to  an  apparent  resistivity  using 
cq  6. 
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Figure  3.  VLF apparent  resistivity  contours  at  an  attitude  of  75  ni  above  a conductive  (a)  and  a 
resistive  (b)  circular  anomaly,  the  anomaly  radius  is  / 50  m and  the  incident  wave  is  from  the  left. 
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RESULTS 

Numerical  integrations 

Equation  21  was  numerically  intr  I to  produce 
apparent  resistivity  contours  at  a standard  survey  alti- 
tude of  75  m above  simulated  circular  anomalies.  Each 
anomaly  was  300  m in  diameter  and  the  resistivity  was 
greater  by  a factor  of  10  than  the  background  resistivity. 
The  areas  of  the  anomalies  were  broken  into  square 
elements  equal  to  100  m2.  This  size  maintained  a dis- 
crepancy of  less  than  1 .76%  between  the  radial  dis- 
tances from  any  two  points  within  the  element  to  the 
airborne  point  of  computation.  The  calculations  used 
the  values  0 = 87°  (Watt  1967),  frequency  = 20  kll/, 
and  assumed  that  the  wave  was  incident  from  the  left, 
as  was  indicated  in  figure  2. 

The  simulated  contours  are  presented  in  I igure  3 for 
a conductive  (I  ig.  3a)  and  a resistive  (Fig.  3b)  anomaly. 

I he  shaded  portions  refer  to  the  anomalous  surface 
areas.  I rom  the  values  above  the  center  of  the  anoma- 
lies, it  is  apparent  how  influential  background  resistivity 
and  altitude  are  upon  the  anomalous  values.  Over  the 
conductive  anomaly,  the  value  ol  660  ohm-m  repre- 
sents a 34%  decrease  from  the  apparent  resistivity  of 
the  background,  yet  it  is  still  6.6  times  larger  than  the 
actual  surface  value.  Over  the  resistive  anomaly,  the 
value  of  260  ohm-m  represents  a 160%  increase  over 
the  background,  yet  is  only  26%  ol  the  actual  surface 


value.  At  an  average  flight  line  spacing  ol  160  m 
(present  practice),  an  aircralt  should  intercept  this 
zone  at  least  within  the  700  ohm-m  contour  of 
I igure  3a  and  approximately  within  the  230  ohm-m 
contour  of  figure  3b. 

I he  slight  asymmetry  in  the  x direction  results 
from  the  inclusion  ol  the  term 

exp  ('/A'q.v'  sin 0) 

lor  the  / ”a  expression.  There  is  no  reason  for  assuming 
that  this  simple  plane  wave  phase  dependency  maintains 
itsell  over  such  a small  anomaly,  but  the  effect  is  only 
marginal  for  Ihe  anomaly  size  considered  in  this  analy- 
sis. 

In  I igurcs  4a  and  4h  are  plotted  curves  that  permit 
the  effects  of  resistivity  contrast  and  area  size  to  be 
separated  for  the  standard  survey  altitude  ol  75  m 
(again,  f- 20  kHz,  0 : 87°,  and  <p  ' 45°  lor  the  entire 
ground  surface).  In  both  I igurcs,  the  background  re 
sistivity  is  given  as  p, , the  anomalous  resistivity  as 
p2,  and  the  independent  vaiiablesarc  anomaly  radius 
ra,  and  resistivity  contrast  (I  p2/p,  Several  compu- 
tations revealed  that  these  curves  could  be  scaled  to 
within  1%  lor  any  resistivity  contrast,  thereby  allow 
ing  this  type  ol  presentation. 

I igure  4a  plots  the  apparent  resistivity  p in  per- 
centage ol  p, , above  the  center  ol  a umdmtive  anomaly 
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Figure  5.  VLF apparent  resistivity  contours  in  the  Tuktoyaktuk  region  of  the  Canadian  Northwest  Territories 
above  a thawed  lake  (a)  and  permafrost  (h).  Map  courtesy  of  Geological  Survey  of  Canada.  Contour  values 
are  in  ohm-m;  shaded  portions  are  areas  of  deep  thaw.  (Aerial  view  of  portion  of  area  h is  shown  on  cover.) 


(p,  > p2).  I inure  4b  plots  pa  in  percentage  of  p2  above 
the  center  of  a resistive  anomaly  (p,  < p2).  Both  fig- 
ures reveal  how  difficult  it  is  to  alter  the  airborne  read- 
ings from  the  background  surface  value  for  an  anomaly 
less  than  1 km  in  diameter.  The  decrease  of  only  50% 
at  a 75 -m  altitude  revealed  in  f igurc  4a  requires  a radius 
of  600  m and  a conductive  contrast  of  10/1.  In  figure 
4b,  a 1 0/1  increase  in  surface  resistivity  for  a radius  of 
600  m results  in  only  a 38%  increase  at  an  altitude  of 
75  m.  Both  figures  imply  that  beyond  the  600-m 
radius  very  little  change  will  take  place  per  100-m  in- 
crease in  radius. 

Any  further  attempts  to  increase  the  radius  beyond 
the  values  given  must  be  limited  by  two  factors.  The 
first  is  that,  as  the  diameter  increases  towards  an 
appreciable  traction  of  a frce-spacc  wavelength  (I  5 km 


at  20  kHz),  a model  that  considers  changes  in  all  the 
field  components  must  be  used.  The  second  is  that, 
in  order  to  completely  alter  the  existing  background 
plane  wave  resistivity  value,  the  integration  must  be 
extended  over  an  infinite  radius. 

Example  no.  I 

I igure  5 compares  the  VI  I airborne  resistivity  con- 
tours over  two  ditferent  sections  of  the  Tuktoyaktuk 
region  of  the  Canadian  Northwest  Territories. 

Ihe  area  has  little  relief  and  alternates  between 
seasonally  thawed  lakes  over  deeper  thaw  /ones  (I  ig. 

5a)  and  permafrost  (I  ig.  5b).  The  airborne  differenti- 
ation of  the  thawed  /ones  is  apparent.  I he  thawed 
zones  in  I igurc  5a  have  an  average  apparent  resistivity 
ol  about  250  ohm-m  interspersed  with  some  depressions 
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Figure  6.  Airborne  and  ground  readings  of  VLF  (18.6  kHz)  apparent  resistivity  (b) 
across  a section  of  the  Goldstream  Valley  (a)  near  Fairbanks,  Alaska.  (After 
Hoekstra  et  at.  1975.)  The  ground  readings  have  been  modified  according  to  eg  7. 
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below  50  ohm-m  and  one  maximum  above  600  ohm-m. 
The  permafrost  (Fig.  5b)  has  an  average  of  about  2000- 
3000  ohm-m  with  several  anomalies  above  5000  ohm-m 
and  some  depressions  near  1 000  ohm-m.  One  of  these  de- 
pressions obviously  correlates  with  the  small  thaw  /one 
shaded  at  the  lower  right  of  Figure  5b  and  will  be  used 
as  an  example,  since  its  dimensions  can  be  verified. 

This  small  thaw  zone  is  represented  by  a closed  1000 
ohm-m  contour  contrasted  against  an  average  value  of 
about  2000  ohm-m.  This  50%  change  in  value  is  found 
at  a radius  of  600  m on  the  0 = 0. 1 curve  of  Figure  4a. 
This  value  of  (3  implies  that  the  actual  resistivity  of  the 
anomaly  is  about  200  ohm-m  and  the  approximate 
radius  is  600,  which  gives  an  approximate  area  of  1.13 
km2 . These  values  compare  favorably  with  the  follow- 
ing observations.  From  I igure  5a  the  predominant 


value  for  a large  thaw  zone  is  250  ohm-m.  From  the 
distance  scale  supplied,  the  small  thaw  zone  in  question 
is  approximately  1 .0  km2  in  area. 

Example  no.  2 

A second  example  is  taken  from  the  work  of 
Hoekstra  (Hoekstra  et  al.  1975)  in  the  Goldstream 
Valley  near  Fairbanks,  Alaska,  shown  in  cross  sec- 
tion in  Figure  6a.  Figure  6b  shows  the  correlation 
between  airborne  and  ground  computations  of  apparent 
resistivity  taken  along  a cross  section  ol  the  valley.  The 
ground  readings  were  modified  using  eq  7.  The  detail 
ol  the  ground  readings  is  lost  in  the  airborne  readings, 
yet  both  sets  of  data  maintain  the  same  approximate 
mean  level  between  800  and  1000  ohm-m.  The  air- 
borne curve  is  extrapolated  between  values  recorded 
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alonp  several  flight  lines  that  traversed  the  length  of  the 
valley  floor  and  nearly  perpendicular  to  the  cross  sec- 
tion given. 

The  ground  data  indicate  that  this  central  portion 
of  the  Goldstream  Valley  was  extremely  resistive  with 
quadrature  ground  values  possibly  exceeding  3000 
ohm-m  against  a background  of  about  1000  ohm-m. 

The  width  of  this  portion  is  about  800  m.  A circular 
approximation  | contours  presented  in  Hoekstra’s 
(Hoekstra  et  al.  1975)  article  indicated  that  high  re- 
sistivities are  not  continuous  along  the  valley  floor)  of 
400-m  radius  with  a resistivity  contrast  ol  3 is  presented 
in  f igurc  4b.  The  figure  shows  that  only  about  57% 
of  the  true  ground  value  is  obtained  in  the  air.  From 
Hoekstra  s airborne  data,  a maximum  value  of  only 
about  1400  ohm-m  is  attained  by  the  airborne  system 
over  this  region,  which  is  47%  of  the  estimated  true 
level. 


CONCLUSIONS 

figures  4a  and  t>  present  a simple  means  of  recover  - 
ing  ground  detail  that  is  lost  by  altitude  in  an  airborne 
VLI  survey.  They  also  permit  preliminary  evaluations 
to  be  made  ol  the  success  of  an  airborne  survey  in  find- 
ing conductive  or  resistive  /ones  covering  less  than  a few 
square  kilometers.  Such  evaluations  would  require  only 
the  desired  value  for  site  resistivity  and  estimates  of  the 
average  values  for  the  region  considered.  Such  estimates 
can  now  be  made  rapidly  with  commercially  available 
and  portable  VLf  field  equipment. 
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